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Abstract 

We present and analyze the simple analytically solvable model of l/f noise, which 
can be relevant for the understanding of the origin, main properties and parameter 
dependencies of the flicker noise. In the model, the currents or signals represented 
as sequences of the random pulses, which recurrence time intervals between transit 
times of pulses are uncorrelated with the shape of the pulse, are analyzed. It is shown 
that for the pulses of fixed area with random duration, distributed uniformly in a 
wide interval, l/f behavior of the power spectrum of the signal or current in wide 
range of frequency may be obtained. 
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1 Introduction 

The origin and omnipresence of 1// noise is one of the oldest problems of the contemporary 
physics. Since the first observation of the flicker noise in the currents of electron tubes by 
Johnson [1] , fluctuations of signals and physical variables exhibiting behavior characterized 
by a power spectral density diverging at low frequencies like l/f have been observed in 
a wide variety of systems [2]. The widespread occurrence suggest that some underlying 
mechanism might exist. However, a fully satisfactory explanation has not yet been found 
and the general theory of 1// noise is still an open question. 

A simple procedures of integration or differentiation of the convenient (white noise, 
Brownian motion or so) fluctuating signals do not yield in the signal exhibiting l/f noise. 
There are no simple linear, even stochastic, differential equations generating signals with 
l/f noise. Therefore, l/f noise is often modeled as the superposition of Lorentzian spectra 
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with a wide range of relaxation times [3]. Summation or integration of the Lorentzians 
with the appropriate weights may yield 1// noise [4]. 

In many cases the physical processes can be represented by a sequence of random pulses. 
Recently, considering signals and currents as consisting of pulses we have shown [5-7] that 
the intrinsic origin of 1// noise may be a Brownian motion of the interevent time of the 
signal pulses, similar to the Brownian fluctuations of the signal amplitude, resulting in 
1//^ noise. 

The model, proposed in [5,6], can be extended taking into account finite duration of 
the pulse. The spectrum of the signal, consisting of the pulse sequences which belong to 
the class of Markov process, was investigated in [8,9]. 

In this article we present a different model of pulses. We consider a signal consisting of a 
sequence of uncorrelated pulses. The shape of the pulses is determined by only one random 
parameter — pulse duration. We will show that by suitably choosing of the distribution 
of the pulse duration the 1// noise can be obtained. 

2 Signal as sequence of pulses 

We will investigate a signal consisting from a sequence of pulses. We assume that: 

(i) the pulse sequences are stationary and ergodic; 

(ii) interevent times and the shapes of different pulses are independent. 
The general form of such signal can be written as 

I{t) = Y,A,{t-t,) (1) 

k 

where functions Ak{t) determine the shape of individual pulses and time moments tk de- 
termine when a pulse occurs. The power spectrum is given by the equation 



S(f) = lim ( ^ 



Jti 



(2) 



where T = tf — ti. Substituting Eq. ([T]) into Eq. ([2]) we have 



Siu) = hm ( - Y: e^-(*^-*^') / ' du f " du'Ak{u)Akiuy'^^^-^'^ ) . (3) 

We assume that functions Ak{u) decrease sufficiently fast when \u\ oo. Since T oo, 
the bounds of the integration in Eq. can be changed to ±oo . We also assume that time 
moments t^. are not correlated with the shape of the pulse A^. Then the power spectrum 
is 

O , , / r+oo r+oo \ 

Siiu) = hm - (e'-^'^-'^'A ( / du ^iw'Afc(w)A,,K)e*"("-" ^ 

^ kk' \J-oo J-OO I 
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After introducing the functions 

/ /•+00 . r+oo . \ 

= ( / rfMAfc(n)e^-" / du'Akiu')e-''^^ ) (4) 

\J—oo J—oo I 

and 

XkA^) = (e-^*'^-*^')) (5) 

the spectrum can be written as 

2 

S{u) = \im -Y,Xk,k'{uj)'^k,k'{uj). (6) 
2.1 Stationary process 

Equation ([6]) can be further simphfied assuming that the process is stationary. In the 
stationary case all averages can depend only on — k'. Then 

'^k,k'H = ^fc-fc'(^)- (7) 

and 

= Xfc-fe'(^)- (8) 

Equation then reads 

2 

^(^) = lim -J2Xk-k'{uj)'^k-k'{uj). 

k,k' 

Changing the variables into k = k' and q = k — k' and changing the order of summation 
we obtain 

2 ^max ^min ^max Q 
g=l fc — ^min 

+^1™^ E E x.M^^M 

Q—^min ^max ^ — ^min Q 
2 ^max 

+ lim — V ^o(t^)- 

/C "'mill 

Introducing = kma.x — ^min we have 

S{u) = 2M^)^ + lim 4 E (z> - I) Rex,(a;)vi>,(cu) (9) 

where 

T^oo \ T 
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i?=lim^— (10) 



is the mean number of pulses per unit time. 

When the sum Ylq=iQ^^Xq{^)'^q{^) converges and T ^ oo then the second term in 
the sum ([9]) vanishes and the spectrum is 

oo 

S{uj) = 2z7^oM+4z7^Rex,M%H (11) 

9=1 

oo 

= 2u X.M^.H- (12) 

q= — oo 

2.2 Fixed shape pulses 

When the shape of the pulses is fixed (/c-independent) then the function '^k,k'{^) does not 
depend on k and k' and, therefore, '^k,k'{^) = ^o,o(^^)- Then equation yields the power 
spectrum 

2 ^ 

S{uj) = ^^0,0(^^)11111 — 2^Xfc,fc'(^) = '^ofl{^^)Ss{uj). (13) 

k,k' 

This is the spectrum of one pulse multiplied by the spectrum of the sequence of 5-shaped 
pulses Ss{uj). It has been shown [5,6] that the spectrum of such a sequence can exhibit 
l//-like behaviour in a broad frequency range if the interevent times Tk = tk — tfc-i follow 
an autoregressive process. 

2.3 Uncorrelated pulses 

When the pulses are uncorrelated and k ^ k' then 

= \{Fk{uj))\' 

where 

/+00 
Ak{u)e'''^du. (14) 
-00 

is the Fourier transform of the pulse A^. When k = k' then 

^o{uj) = {\Fk{uj)\') . 
/^From Eq. (fTTj) we obtain the spectrum 

00 

Siu)=2u(\Fkiu;)\') + Au\{Fk{cu))\'J2^eXqiu;). (15) 

q=l 

When the interevent times Tk = tk — tfc-i are random and uncorrelated then 

X.H = (e^-(*'=+.-*^)) = (e^-^y = xAujy. (16) 
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^From Eq. (|T^ we obtain 



2z7(|FfcHn+4z7|(Ffc(u;))|2Re 



1 - xA^)' 



(17) 



Here 









a;=0- 



If the occurrence times of the pulses tk are distributed according to Poisson process then 
the interevent time probabihty distribution is '^{t) = ^e~^. The characteristic function 
obeys the equahty Re i^^^*^^) = and the spectrum is 

S{u) = 2u{\Fk{uj)\^). (19) 
We will investigate this case more deeply. 



3 Pulses of variable duration 

Let the only random parameter of the pulse is the duration. We take the form of the pulse 

as 

A,{t) = T,^A (^) , (20) 

where Tk is the characteristic duration of the pulse. The value P = corresponds to fixed 
height pulses; /? = — 1 corresponds to constant area pulses. Differentiating the fixed area 
pulses we obtain f3 = —2. The Fourier transform of the pulse fl20l) is 

T^A — ) e'^'dt = +^ / A(M)e^-^'="rfM = +^F(u;Tfc). 

-oo \1 k / J —oo 

iFTom Eq. ffTI?|) the power spectrum is 

S{uj)=2u{T',^^'\F{u;n)\'). (21) 
Introducing the probability density P(Tk) of the pulses durations Tk we can write 

S{u) = 2v r Tlf'^''\F{ujTk)\^P(Tk)dTk. (22) 

JO 

If PiTk) is a power-law distribution, then the expressions for the spectrum are similar for 
all (5. 
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3.1 Spectrum at small frequencies uj 

For small frequencies we expand the Fourier transform of the pulse into Taylor series. The 
first coefficients are 

= "ij^^'""^^^' = ^^^^ 

where 

/+00 
A{t)dt (24) 
-oo 

is the area of the pulse, 

1 f + OO 1 f + OO 

{t) = - tA{t)dt, (t^) = - / t^A{t)dt. (25) 

Qj J— oo Qj J —oo 

Then the spectrum from Eq. ([221) 

roc 

S{lu) ^ 2va^ / Tf +2(1 - At^u;'T^)P{Tk)dn , 
Jo 

where At^ = (t^) - (t)^ We obtain 

Siu) = 2ua'{T'/^'){l - At'u'iT'/^')), (26) 

where 



{n)= nP{T,)dTu. (27) 







3.2 Power-law distribution 

We take the power-law distribution of pulse durations 



P{Tk) = \ S^^^' ' (28) 

1 0, othervise. 



IFtoiu. Eq. (|22|) we have the spectrum 

a + 1 



max min ^ 

2u{a + 1) /■'^^. 



UJ' 



a+2/3+3|'7^a+l _ 7^"+l^ /, ,^ 



When a > — 1 and t^t^ ^ cx; ^ tJ— then the expression for the spectrum can be approxi- 

^ max -t min 

mated as 

I-' max -'min J •^'^ 
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If a + 2/3 + 2 = then in the frequency domain tJ— <^ a; <^ tJ— the spectrum is 



max 



^ ^7^^^^^ L \Fiu)\-'du. (30) 

\ max ir 



' mm 



We obtained 1/f spectrum. The condition a + 2/5 + 2 = is satisfied, e.g., for the fixed 
area pulses (/3 = —1) and uniform distribution of pulse durations or for fixed height pulses 
(/3 = 0) and uniform distribution of inverse durations 7 = TjT'^ , i.e. for P{Tk) oc T^^ . 
If a + 2/3 + 4 = then in the frequency domain tJ— <^uj <^ tJ— the spectrum is 

J- max J- mill 

I-' max -'min J "^'^ " 

Such a spectrum can be obtained after differentiation of the signal exhibiting 1/f spectrum. 



4 Example 

4.1 Rectangular pulses 

As an example we will obtain the spectrum of rectangular constant area pulses. The 
duration of the pulse is T^. The Fourier transform of the pulse is 



1 _ pii^Ti 



_1_ ^^2sin(^f^) 



F{uTk) = a / due''^^'^^ = a . ^ = ae'— (32) 



Then the spectrum according to Eqs. fl22l) . flzHjl and fl32|) is 

4z7a^(a + l)(r^-i-r^rn^) , 4z7a2(a + l) 



V" -^/V-^max -^mm / ^ V-^max mm / 

X Re{z^'" (r(a - l,2cuT^ax) - r(a - l.iujT^^))] , (33) 

where r(a, 2;) is the incomplete gamma function, r(a, 2;) = u"'~^e~^du. 

When — 1 < a < 1 then the term with T{a — l,iu!T^ax) is small and can be neglected. 

We also assume that Tmin <^ T^s^x and when a < — 1 we neglect the term I tj™- . Then 

\ max / 

we have 

^, . 4z/a^(a + l) /TT .xAt./ .X /O.N 

For a = we have the uniform distribution of the pulses duration. Using the result of the 
limit 

limcos (^^(a- 1)^ r(a- 1) = (35) 

we obtain 1/f spectrum 

Siu) ^ (36) 

max 

The spectrum was also obtained from numerical calculations. Typical signal for rect- 
angular fixed area pulses is shown in Fig. [T] and the power spectrum in Fig. [21 
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Figure 1: Typical signal consisting from the fixed area rectangular pulses with uniformly 
distributed durations. The time intervals between the pulses are distributed according to 
Poisson process with the average f = 5. The used parameters are T^in = 0.01, Tmax = 100 
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Figure 2: The spectrum of the signal consisting from the fixed area pulses with uniformly 
distributed durations. The used parameters are the same as in Fig. [TJ The dashed line 
corresponds to the spectrum obtained according Eq. flHUj) 
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Figure 3: The spectrum of differentiated signal consisting from the fixed area pulses. The 
used parameters are the same as in Fig. [TJ The dashed line corresponds to the spectrum 
obtained according Eq. ( HOl) . 



4.2 Differentiated rectangular pulses 

After differentiation of the signal the shape of the pulse becomes 



The Fourier transform of the pulse is 

F(cjTfc) = a(l - e*^^'= 

Then the spectrum is 



—2iae^ 2 gin 



ujTi, 



(37) 



(3J 



V" y \ max mm / ^ \ max mm / 

X Re {i^~°'{T{a - 1, iujT^^^) - T{a - 1, icuTmin))) . 



When a = we obtain 



277 z/g^ 

max 



(39) 
(40) 



The spectrum, obtained numerically, is shown in Fig. [31 
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5 Conclusion 



We investigate signals consisting of a sequence of uncorrelated pulses with random du- 
rations. By suitably choosing distribution of the pulse duration the 1// noise can be 
obtained. Signal of fixed area pulses yields 1// noise when width of the pulse is uniformly 
distributed in a wide interval. The spectrum is given by Eq. (IHUI) . This conclusion does 
not depend on particular shape of the pulse. For the fixed amplitude pulses 1// spec- 
trum yields when the inverse duration of the pulses 7 = Tj^^ is distributed uniformly or 
PiTu) oc T,-2 
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